one  hiLE  copy 


NPS012-87-006 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


C 

ELECTE 

MAY  1  Ot988 


THE  DIFFUSION  OF  VORTICITY 
FROM  A  PLANE  BOUNDARY 
by 

A.  V.  Hershey 
September  1987 


Approved  for  public  release;  distribution  i 

Prepared  for:  Naval  Postgraduate  School 
Monterey,  CA  93943 

s  unlimited 

88  S 

09 

NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  CA  93943 


RADM  R.C.  AUSTIN  K.T.  MARSHALL 

Superintendent  Acting  Provost 


This  report  is  a  contribution  to  numerical  ship  hydrodynamics. 
Prepared  by: 


K.\y.  Ha 

A.V.  HERSHEY  \ 


Research  Affiliate 


Reviewed  by: 


D.G.  WILLIAMS 
Director,  Computer  Center 


Ij  REPOHI  StCURllY  LlASSihCAIlON 

Unclassified 


i*  SECURITY  CLA^SihLATlUN  AUIhORIIY 


20  OECt-ASSiTiCATlON  /  UOWNCRAOING  SCHEOULk 


4  PkKTOHIVliNu  UKoA,Si2AllON  HkPOKI  NUIVItlkK(S) 

NPS012-87-006 


REPORT  DOCUMENTATION  PAGE 


ID  KESIHikllVk  MAhii.iNkiS 


J  OlSTMlBuTIUN/ AVAiLAdlUlY  OP  REPORT 

Approved  for  public  release;  distribution  is 
uni imited. 


)  MUNilOHiNG  ORoANi^ATiON  RkPOKT  NuMdkR(S) 


Dj  NAME  OP  PERPOKMiNO  ORUANiZAtlON  bO  OPPlCE  SYMBOL  /p  NAME  OP  MONITORING  ORGANISATION 
S  (If  ippinibit) 

Naval  Postgraduate  School 


6<  ADDRESS  (City.  StRte.  jnd  /iPCodf) 

Monterey,  CA  93943 


a<t  NAME  OP  PuNOiNG/SPONSORiNG 
ORGAMSATiON 


B<  ADDRESS  (City.  Stite,  dnd  iIPCodt) 


Z  I  )i  TiTlE  tindude  Security  Clitsihcstion) 

u  I 


/b  ADDRESS  (City,  Sure,  jnoSlPComi 


80  OPPlCE  SYMBOL  |  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  ipphtJbl*)  I 


!0  sottRCE  OP  Funding  numbers 


PROGRAM 


PROJECT 

task 

NO 

NO 

AORK  UMT 
ACCESSION  NO 


The  Diffusion  of  Vorticity  from  a  Plane  Boundary 


IS  personal  AuTmOR(S) 

A.V.  Hershe 


134  type  OP  REPORT 

Final 


16  Supplementary  notation 


I3b  TIME  COVERED 
FROM  12/86  TO 


14  DATE  Of  REPORT  (Y*4r  Month.  04y3  15  PAGE  COUNT 

September  1987  24 


COSATi  CODES 


GROUP 


18  SUBlECT  TERMS  (Continue  on  reverie  it  neteiiery  end  Rtentity  by  b/ock  number) 

SUB  GROUP  Fourier  Analysis;  Hydrodynamics  ,  ^5=*^' 

Numerical  Analysis  - 


19  abstract  (Continue  on  reverie  if  neceiiery  end  Kttnvfy  by  bttxfc  number) 


Any  discontinuity  in  flow  at  a  plane  boundary  is  dispersed  by  the  diffusion  of  vorticity 
into  a  boundary  layer  and  v/ake.  Experimental  measurements  in  the  past  have  been  correlated 
with  empirical  formulations.  A  generalization  gives  both  the  horizontal  component  and  the 
vertical  component  of  the  mean  velocity.  In  the  laminar  sublayer  the  velocity  is  a  solution 
of  the  diffusion  equation.  In  the  turbulent  boundary  layer  the  velocity  can  be  expressed  by 
a  Fourier"  integral.  In  the  free  stream  there  is  a  vertical  persisteVice  of  velocity.  The 
computation  of  velocity  is  provided  by  subroutines.  /,  / 


20  Distribution  I  AVAILABILITY  of  abstract 
□  uNClASSifiEDiiiNLiMiTED  □  SAME  AS  RPT  □  DTiC  USERS 


224  NAME  Of  RESPONSIBLE  INDIVIDUAL 

A.V.  Hershe 


21  ABSTRACT  security  CLASSifiCATlON 

Unclassified 


22b  telephone  (mcAidt  *n»  Cod*)  22c  Of  fKE  SYMBOL 


iftSiMSiStFbl.' 


DDFORM  1473, 84  MAR 

83  APR  edition  m4y  be  uted  until  eiheuitcd 

All  other  editioni  ere  obtolett 

NPS012-87-006 


September  1987 


THE  DIFFUSION  OF  VORTICITY 


FROM  A  PLANE  BOUNDARY 


By 

A.  V.  HEKSHEY 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY.  CALIFORNIA  93943 


Aooesslon  For 


HTIS  fflKUbI 
DIIC  TAB 
Unannounoed 
Justlfioatloa. 


By_ 


Distribution/ 


Availability  Codes 
jAval^  and/or 
Special 


Dlst 


m 


Approved  for  public  release,  distribution  unlimited 


TABLE  OF  CONTENTS 


Page 

TABLE  OF  CONTENTS  .  i 

ABSTRACT  .  ii 

INTRODUCTION  .  1 

EQUATION  OF  CONTINUITY .  1 

EQUATION  OF  MOTION .  2 

BLASIUS  PROFILE .  3 

TURBULENT  PROFILE .  5 

LINE  FLUX  PROFILE .  7 

NONSLIP  BOUNDARY .  11 

FOURIER  TRANSFORM .  12 

INCOMPRESSIBLE  FLUID  .  15 

DIFFUSION .  17 

DISCUSSION .  18 

CONCLUSION .  19 

BIBLIOGRAPHY  .  20 

DISTRIBUTION .  21 


ABSTRACT 


Any  discontinuity  in  flow  at  a  plane  boundary  is  dispersed  by  the  diffusion  of  vorticity 
into  a  boundary  layer  and  wake.  Experimental  measurements  in  the  past  have  been 
correlated  with  empirical  formulations.  A  generalization  gives  both  the  horizontal 
component  and  the  vertical  component  of  the  mean  velocity.  In  the  laminar  sublayer 
the  velocity  is  a  solution  of  the  diffusion  equation  In  the  turbulent  boundary  layer 
the  velocity  can  be  expressed  by  a  Fourier  integral.  In  the  free  stream  there  is  a 
vertical  persistence  of  velocity  The  computation  of  velocity  is  provided  by  subroutines 


INTRODUCTION 


The  drag  on  a  ship  is  partly  a  viscous  resistance  and  partly  a  wave  resistance.  The 
viscosity  of  the  fluid  causes  a  transport  of  momentum  from  the  ship  to  the  fluid  The 
momentum  is  located  in  a  boundary  layer  at  the  surface  of  the  ship,  and  in  the  wake 
behind  the  ship  The  boundary  layer  is  laminar  at  the  bow,  but  the  boundary  layer  is 
turbulent  at  the  stern  For  a  nonslip  boundary  condition  the  velocity  relative  to  the 
ship  IS  zero  at  the  surface  and  is  the  velocity  for  free  flow  far  from  the  surface.  A 
knowledge  of  the  velocity  may  be  found  from  an  analysis  of  the  analogous  problem 
in  the  flow  over  a  plane. 

When  a  thin  plate  moves  edgewise  through  a  viscous  fluid  the  plate  entrains  fluid 
if  there  is  a  nonslip  boundary  condition  at  the  surface  of  the  plate.  Each  differential 
element  of  surface  of  the  plate  sets  up  a  current  which  trails  downstream  and  spreads 
out  by  diffusion.  The  velocity  at  the  plate  is  the  accumulated  sum  of  velocities  in 
currents  which  have  been  created  upstream. 

The  boundary  layer  of  a  plate  has  been  the  subject  of  many  investigations’"’  There 
are  three  principal  methods  of  analysis  The  first  method  is  statistical.  It  leads  to  a 
system  in  which  the  number  of  variables  is  greater  than  the  number  of  equations. 
Various  schemes  for  closure  have  been  proposed  The  second  method  is  polynomial 
The  velocity  is  expressed  as  a  power  polynomial  in  a  limited  region  of  a  laminar 
boundary  layer.  The  third  method  is  spectral  The  velocity  is  expressed  as  a  Fourier 
integral  The  rate  of  change  of  the  Fourier  amplitude  is  determined  by  an 
integro-differential  equation  The  evolution  of  a  velocity  distribution  from  an  initial 
distribution  is  determined  uniquely  by  the  integro-differential  equation. 

The  Fourier  analysis  has  been  applied  to  homogeneous  isotropic  turbulence  Contact 
of  the  fluid  with  a  solid  boundary  is  the  origin  of  real  turbulence. 

A  comprehensive  source  of  information  about  boundary  layers  is  Boundary -Layer 
Theory  by  Schlichting’  The  Fourier  analysis  of  flow  in  a  fluid  has  been  investigated 
by  Batchelor®  and  by  Orszag®. 


EQUATION  OF  CONTINUITY 

A  fluid  consists  of  particles  with  random  velocities.  The  number  of  particles  per 
unit  volume  determines  a  median  density  p.  Across  any  mathematical  boundary  there 
IS  a  flux  of  particles  from  both  sides.  A  boundary  which  is  moving  at  such  a  speed 
that  the  flux  is  zero  defines  a  median  velocity  v. 

In  a  continuous  fluid  the  law  of  conservation  of  mass  requires  that  within  any 
mathematical  boundary  the  mass  density  p  and  the  particle  velocity  v  are  related  by 
the  equation 


f  ^ 

J  dt 


dr  -* 


j  pv  ds  =  0 


(1) 


where  t  is  the  time,  dr  is  a  volume  element  within  the  boundary,  and  ds  is  a  surface 
element  on  the  boundary.  Application  of  the  Gauss  theorem  leads  to  the  equation 


dt 


+  V(pv)  = 


dt 


dp 

+  v  Vp  +  pV  v  -•  —  +  pV-v  =  0 
dt 


(2) 


The  divergence  of  velocity  is  the  rate  of  expansion  of  the  fluid 


Sn: 


EQUATION  OF  MOTION 

The  bulk  modulus  a,  is  the  rate  of  decrease  of  pressure  per  unit  rale  of  expansion. 
A  pressure  pulse  is  propagated  through  the  fluid  with  the  speed  of  sound  (v  p)‘ 

In  a  continuous  fluid  the  gradient  of  velocity  satisfies  the  identity 


Cv  =  |(Vv  -  rv)  H-  |(Uv  +  Tv) 


(3) 


where  V'v  is  the  transpose  of  Vv  The  antisymmetric  part  is  the  rate  of  rotation  fJ  and 
the  symmetric  part  is  the  rate  of  strain  6.  There  would  be  spinup  of  fluid  particles 
if  stress  were  not  symmetric.  In  an  isotropic  fluid  the  pressure  is  isotropic  and  stiear 
stress  IS  proportional  to  shear  rate  The  strain  rate  is  partly  an  isotropic  expansion 
rate,  which  is  the  trace  of  the  strain  rate,  and  is  partly  an  anisotropic  distortion  rate, 
which  IS  the  shear  rate  The  shear  stress  arises  from  the  anisotropic  distortion  rate. 
The  stress  E  is  given  by  the  equation 


E  =  -  p  I  -  r-v  I  ^  2p  e 

where  p  is  the  pressure,  and  p  is  the  viscosity. 

The  force  on  any  mathematical  boundary  is  given  by  the  Gauss  theorem 


J  T  E  dr  =  E  ds 


(4) 


(5) 


where  dr  is  a  volume  element  within  the  boundary  and  ds  is  a  surface  element  at  the 
boundary  The  equation  of  motion  per  unit  volume  is  given  by  the  equation 


dv 

p  —  ^  VS 
dt 

As  a  consequence  of  the  identity 

V-V'v  =  VV'V 

the  equation  of  motion  is  the  Navier-Stokes  equation 

dv 
dt 

The  kinematic  viscosity  v  is  given  by  the  equation 

M 

1/  =  — 


vVv  = 


,  f  ^  1  M 

J  P  ^  P 


VV-v  *  -  V  Vv 
P 


(6) 


(V) 


(8) 


(9) 


For  hydrodynamics  the  kinematic  viscosity  is  constant  Differentiation  throughout  the 
equation  of  motion  leads  to  the  equation 


a 


v  V(V-v)  t  (Vv-Vv)  =■-  -  V  V 


J  p  ^ 


3  VV(V-v) 
P 


(10) 


where  a  tensor  enclosed  in  parentheses  is  the  contraction  of  the  tensor  by  internal 
scalar  multiplication  This  equation  is  a  Poisson  equation  for  the  determination  of 
pressure  Differentiation  througtioul  the  equation  of  motion  leads  to  the  equation 


dt 


(V- v)  *  v-V(V»v) 


jVvVv]  =  -  V  V(V- v) 
P 


(11) 


wtiere  a  tensor  enclosed  in  brackets  is  the  contraction  of  the  tensor  by  internal  vector 
mulliplic alion  This  equation  is  a  diffusion  equation  for  the  determination  of  vorticity. 


2 


RLASirS  PROFILE 


For  two-dimensional  flow  it  is  fashionable  to  define  Cartesian  coordinates  x.y.z 
such  that  X  is  parallel  to  the  plane  in  the  direction  of  flow,  y  is  perpendicular  to  the 
plane,  and  c  is  parallel  to  the  plane  in  the  direction  perpendicular  to  the  flow  Let 
i,j.k  be  unit  vectors  in  the  directions  of  increasing  x.y.z  The  particle  velocity  v  is 
given  by  the  equation 

V  =  ui  -  i  j  (12) 


where  the  components  u.  v  are  zero  at  the  plane  but  approach  U.  0  with  increasing  y. 
Continuity  in  an  incompressible  fluid  is  expressed  by  the  equation 


du  dv 
dx  dy 


(13) 


The  continuity  equation  is  satisfied  when  the  components  of  velocity  are  derived  from 
a  stream  function  in  accordance  with  the  equations 


u  = 


d\l/ 

% 


(14) 


This  follows  from  the  principle  that  a  second  order  derivative  is  independent  of  the 
order  of  differentiation 

For  stationary  flow  the  derivative  with  respect  to  t  is  zero  Far  downstream  derivatives 
with  respect  to  x  approach  zero,  and  the  NavK .  .'■’tokes  equation  is  reduced  to  the 
equation 


du  du  d'u 
dx  ^  dy  '  dy^ 


(15) 


Far  downstream  the  profile  of  the  boundary  layer  approaches  a  constant  limit. 

The  limiting  profile  is  given  by  a  relation  between  two  dimensionless  variables 
^  and  which  are  defined  by  the  equations 

(77  V  — 

?  =  \ I - ^  V'=vi'f'x<-  (16) 

\  e  Vx 


Differentiation  leads  to  the  equations 

u=(7<-'  ^  =  (17) 

Then  substitution  m  the  reduced  equation  of  motion  leads  to  the  Blasius  equation 

cr  +  2r'  =  0  (18) 

Boundary  conditions  for  this  differential  equation  are  given  by  the  equations 

<■  =  <^’ =  0  at  ^  0  (19) 

<■'-1  at  (20) 

In  the  original  derivation  by  Blasius'  the  variables  were  so  defined  that  approached 
a  limit  of  2  The  derivation  herein  follows  the  derivation  by  Schlichting’  where  the 
variables  are  so  defined  that  approaches  a  limit  of  unity. 


An  ascending  senes  is  given  by  the  equation 


I  (^^1) 

m  =  0 

This  series  satisfies  the  boundary  conditions  at  ^  =  0.  Substitution  in  the  differential 
equation  leads  to  the  recurrence  equation 

m-1 

Z  (3m  -  3k  -  2)(3m  -  3Jt  -  +  2(3m)(3m  +  l)(3m  +  2)a^  =  0  (22) 

t=o 

The  recurrence  is  started  with  a  value  of  Oq  which  is  adjusted  by  trial  to  make  the 
series  meet  the  boundary  conditions  at  ^ o°  The  computed  value  of  Qq  is  0.16603. 
The  ratio  between  successive  coefficients  of  the  series  tends  to  a  constant  limit  with 
increasing  order.  It  is  possible  to  estimate  a  remainder  after  a  finite  number  of  terms 
on  the  basis  of  the  geometric  series  In  no  case  can  the  series  be  used  for  (  greater 
than  5.69,  where  the  ratio  between  terms  becomes  unity 

For  large  values  of  ^  the  value  of  ^  is  the  sum  of  a  linear  term  and  a  correction 
Correct  to  first  order  the  correction  is  a  solution  of  the  equation 

(i  +  c)r-2r'  =  o  (23) 

where  c  is  an  arbitrary  constant  The  computed  value  of  c  is  -1.72077,  The  solution 
of  the  differential  equation  is  given  by  the  equations 

^  ^  +  c  -  A(^  +  c)  r  d(,  -  Z4 

=  I  -  A  j 

where  A  is  an  arbitrary  constant.  The  arbitrary  constant  is  selected  so  as  to  make 
the  corrected  linear  terms  coincide  with  the  series  expansion  at  ^  =  5.0  The  integral 
in  the  correction  is  given  by  the  equation 

rtCO 

J  [  1  -  erf(|((  +  c))]  (26) 

where  erf  is  the  error  function. 

Previous  analyses  and  computations  on  the  Blasius  profile  are  summarized  by  the 
data  in  a  table  in  a  paper  by  Howarth®,  which  is  the  basis  for  the  Table  7  1  in  the 
text  by  Schlichting’  That  table  of  data  is  reproduced  by  the  following  subroutine 

S'JpROUTiiJE  BLSSPF  (AU,  AN,  AX,  AY,  FU.  FV) 


FORTRAN  SUBROUTINE  FOR  BLASIUS  PROFILE 


(24) 

(25) 


The  free-stream  velocity  U  is  given  in  argument  AU,  and  the  kinematic  viscosity  u 
IS  given  in  argument  Af;  The  coordinates  x.y  are  given  in  the  arguments  Av  a>  Series 
expansions  and  error  integrations  are  used  in  the  evaluation  of  the  Blasius  profile. 
The  components  u,  i'  are  stored  in  functions  FU.  FV. 


TL  RBl’LENT  PROFILE 


Far  downstream  the  flow  in  the  boundary  layer  is  turbulent  There  have  been  many 
experiments  on  the  turbulent  boundary  layer  In  the  absence  of  a  fundamental  theory 
the  experiments  have  been  correlated  with  empirical  formulations.  A  celebrated 
formulation  is  the  logarithmic  law,  which  expresses  mean  velocity  as  the  logarithm 
of  the  distance  from  a  wall  The  logarithm  cannot  be  used  at  the  wall  where  velocity 
IS  zero  and  the  logarithm  is  -«>.  or  at  infinite  distance  from  the  wall  where  velocity 
is  finite  and  the  logarithm  is  +«>.  Corrections  have  been  published  by  Reichardt*  and 
by  Thompson*  The  argument  of  the  logarithm  is  incremented  by  unity  and  the  logarithm 
IS  blended  with  the  free-stream  velocity. 

The  Reynolds  number  is  defined  by  the  equation 


where  V  is  the  free-stream  velocity,  v  is  the  kinematic  viscosity  and  x  is  the  distance 
downstream.  The  drag  D  on  the  wall  is  given  by  the  equation 


D  =  i  pU^ 


cx 


(log  R^)‘ 


(28) 


where  the  constant  c  is  defined  by  the  equation 

c  =  0.455  (log  10)2  *®  (29) 

This  formulation  is  equivalent  to  Equation  21.16  in  the  text  by  Schlichting’.  A  shear 
velocity  u.  is  defined  by  the  equation 


(30) 


where  t  is  the  shear  stress  at  the  wall,  and  p  is  the  density.  The  shear  stress  t  is 
given  by  the  equation 


— 

dx 


(31) 


The  boundary  layer  thickness  <5  is  determined  by  the  equation 


u*6  Ri 

V  (log/?,)2** 


(32) 


This  equation  is  equivalent  to  Table  21.1  in  the  text  by  Schlichting’.  Finally  the  mean 
velocity  u  is  given  by  the  equation 


ru.  / 

u  =  —  log  I 


1  +  K 


u.y 


:i  -7)U 


(33) 


where  x.  is  a  constant  and  y  is  the  blending  function  The  experimental  value  of  k  is 
given  by  the  equation 

1C  =  0  40  (34) 


Various  empirical  schemes  for  the  blending  function  have  been  tried,  but  a  more 
logical  basis  for  the  blending  function  would  be  the  error  function  To  within 


expermienlal  error  the  blending  function  can  be  expressed  by  the  equations 


77  =  TT 


j  erf(  1;  *  5 )  -  erf(77  -  |) 
^  erf(") 


(35) 


Then  the  blending  function  is  in  the  range  1  >  7  >  0  while  y  is  111  the  range  0  <  y  < 
Differentiation  with  respect  to  x  leads  to  the  equation 


ipU‘ 


58 


(log 

and  further  differentiation  leads  to  the  equation 

dr  ,  2  58c 

~  ^Uog 


log 


I  - 


3  58 
log  R^ 


Then  the  derivative  of  u.  is  given  by  the  equation 


du,  .  1  dr 

dx  ^  r  dx 


Differentiation  with  respect  to  x  leads  to  the  equation 


L- 


dd 


dx  u.(logR^) 


3  58 


1  - 


log  Rj.  -  |{2.58)- 


1.58 
log  Rj. 


1 


58 


log  R, 


Then  the  derivative  of  7  is  given  by  the  equation 


^  =  v'n  -  - 
dx  <5 


2  2  ^ 

(5  di 


erf(i) 

Finally  the  derivative  of  u  with  respect  to  x  is  given  by  the  equation 


(36) 


(37) 


(38) 


(39) 


(40) 


du 

dx 


dx 


1/ 


du. 


u.y  dx 


dx 


(41) 


Although  the  derivatives  are  given  by  empirical  functions  with  finite  numbers  of  terms 
their  integration  with  respect  to  y  would  lead  to  an  unlimited  number  of  terms  In  a 
practical  integration  with  respect  to  y  the  integrand  is  best  approximated  by  a  power 
polynomial  in  y  and  integrated  term  by  term  A  discrete  set  of  integrands  is  converted 
into  a  discrete  set  of  coefficients  by  1 1  -point  Lagrange  interpolation  and  the  integration 
IS  completed  coefficient  by  coefficient 

The  stream  function  V  is  given  by  the  equation 

V"  =  r  udy  (42) 

*'o 


6 


KXjnonyrkKWTiLtv 


(43) 


$ 


V 


i 

r^l 
*•  *• 

I 


i 

*.  N 

*»yi 

'’-‘*•3 


’  Sv 


and  the  vertical  coinfwnenl  ?'  is  given  by  the  equation 

1’^  du 

»'  =  -  dy 

Jo 

The  integrations  are  performed  with  the  aid  of  the  following  subroutine. 

subroutine  TBLN'='t  (AU,  an,  ax,  -AY.  FU.  FV) 


FORTRAN  SUBROUTINE  FOR  TURBULENT  PROFILE. 


The  free-stream  velocity  U  is  given  in  argument  AU,  and  the  kinematic  viscosity  u 
is  given  in  argument  AN,  The  coordinates  x.  y  are  given  in  the  arguments  AX.  AY. 
Empirical  formulations  and  Lagrange  interpolation  are  used  in  the  evaluation  of  the 
turbulent  profile  The  components  u.i  are  stored  in  the  functions  Fj. 

LINE  FLUX  PROFILE 

The  usual  assumption  is  made  that  pressure  is  constant  over  a  flat  plate  in  a  steady 
flow  parallel  to  the  plate  Let  the  free-stream  velocity  -^f'i  be  disturbed  by  a  small 
counter  velocitv  v.  The  Fourier  transform  can  be  applied  to  a  function  which  is  zero 
everywhere  except  at  the  origin  Then  Fourier  integration  leads  to  the  equation 


That  this  expression  for  velocity  has  zero  divergence  can  be  verified  directly  by 
differentiation.  Furthermore  it  can  be  derived  from  the  stream  function  in  the 
equation 


^  27t  iac 


At  X  =  0  an  application  of  the  Euler  theorem  and  an  integration  of  the  sine  quotient 
function  stiows  that  V"  ~  ~  |  foi"  all  y  >  0  and  |  for  all  y  <  0  Thus  the  flux  in  the 

counter  flow  is  unity  at  x  =  0  The  curl  of  the  velocity  is  given  by  the  equation 


Vxv  =  + 


(A"  -  x") 


To  within  small  quantities  of  second  order  the  vorticity  satisfies  the  differential  equation 


ax  \ 


d^y  d^y 


Substitution  of  the  Fourier  integral  into  the  differential  equation  leads  to  the  equation 

f'A --  I- (48) 
which  may  be  solved  by  the  quadratic  rule  to  give  the  equation 


2i'  Wui"'  r  ^  \  L'^  * 


Ttie  parameter  A  is  an  even  function  of  the  parameter  k  The  radical  is  negative 
upstream  and  is  pc^sitive  downstream  Ottierwise  ttie  integrals  do  not  converge 


When  y  is  increased  to  infinity  the  integration  witti  respect  to  k  makes  a  sicnifie.int 
contribution  to  ttie  integration  only  where  k  0. 

Upstream  the  parameter  \  is  given  by  the  limit 


V 


X  -  - 

I' 

Thus  the  stream  function  ^  is  given  by  the  limit 

(  j! 

and  the  velocity  v  is  given  by  the  limit 

.  V  -Ui 


J 


(ii  -  0)  (50) 


(y  ^  ±  o=)  (51) 


(y  ^  i  (52) 


The  vertical  component  of  velocity  persists  with  increasing  y.  but  diminishes  rapidly 
with  distance  x  upstream 

Downstream  the  parameter  X  is  given  by  the  limit 


X  -  0 

Thus  the  stream  function  ^  is  given  by  the  limit 


{k  -  0)  (53) 


(y  ->  ±  ~)  (54) 


and  the  flux  in  the  counter  current  is  everywhere  unity.  It  is  independent  of  distance 
X  downstream. 

In  the  integrals  for  the  components  of  velocity  the  integrands  are  the  products  of 
the  monotomc  factors 


2tt 


and  the  oscillatorv  factor 


ZniK 


(55) 


(56) 


The  integration  through  any  number  of  cycles  of  the  oscillatory  factor  can  be  completed 
if  the  monotomc  factors  are  expressed  as  power  series  in  k.  The  range  of  the  power 
series  is  limited  by  the  presence  of  the  radical 


\(^) 

in  the  parameter  X  For  small  values  of  k  the  radical  is  expressed  by  the  equation 

(-l)’"*'(27n)' 


(£) 


n 

=  —  v> 


\\2n/  "  2r  \  U 

This  senes  converges  only  when  k  meets  the  limitation 


2i’ 


U 

Ti' 


(58) 


(59) 


The  senes  is  an  even  function  of  k  For  large  values  of  k  the  radical  is  expressed  by 
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the  equation 


Tins  senes  convert 


\ 


(£)■ 


V 


,*^0  (2m  -  n2^"'(m 

only  when  k  meets  the  limitation 

2i'  1  2i' 

- <  -  <  +  — 

L  K  U 


(60) 


(61) 


The  series  is  an  odd  function  of  1/  a;.  Efficient  evaluation  of  the  series  is  only  possible 
if  the  ranges  of  their  arguments  are  much  less  than  the  limits  of  their  convergence. 
Ttie  monotonic  factors  are  expanded  in  each  of  a  sequence  of  intervals  of  limited 
range 

The  value  of  x  in  the  first  interval  is  given  by  the  equation 


K  =  1]  ~  d 


(62) 


where  ij  is  the  center  of  expansion  and  0  is  the  variable  of  expansion.  The  variable  0 
IS  given  by  the  equation 

0  =  cu  (-  1  g  u  S  1)  (63) 

where  e  is  half  the  range  of  expansion  and  u  is  a  variable  of  interpolation.  The  first 
interval  straddles  the  origin  where  77  =  0.  then  in  subsequent  intervals  rj  is  incremented 
by  2€.  Thus  the  monolonic  factors  are  approximated  by  the  series 


I  ar„u’ 


\ 

K 


m*0 


Required  for  the  computation  of  velocity  are  the  integrals 


c  c 


Jo 


U  £' 


‘  du 


€  C 


n.  ptfyu 


Required  for  the  integration  is  the  recurrence  equation 


■*1  1 

u 

Jo 


du  = 


fu”' 

—  - 

m 

.  Jo 

ify  Jo 

The  recurrence  is  started  with  the  initial  integral  in  the  equation 


/■ 


du 


L  ity  Jo 


The  recurrence  is  cycled  in  ascending  order  if  cy  satisfies  the  inequality 

icyi  g  17 


(64) 


(65) 


(66) 


(67) 


(68) 


Otherwise  the  recurrence  is  cycled  in  descending  order  with  an  initial  approximation 
for'  m  64 

The  value  of  a'  in  the  last  interval  is  given  by  the  equation 


6 

K  = - 

U 


(-  1  g  u  S  4 1)  (69) 


where  r5  is  the  limit  of  integration  The  last  interval  straddles  the  point  where  ^  0 

( Oiitiriuity  through  the  point  is  achieved  by  giving  to  the  radical  in  X  the  same  sign 
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as  the  sign  of  a  The  monolonie  factors  are  approximated  by  ttie  series 


n 

4 

% 


f 

I 

I  1 

wl' 

w!i 

d'St 


•si!: 

ni 

■»! 


•'y 

•' y 
*i»y 

jCV 

iW 

.  .f 

*w 

a 


'(X-<)X  _ 


Then  the  component?  of  velocity  are  expressed  in  terms  of  the  integrals 


£  c„  r(-^)" 

m  =  0  Je  ' 


^,-«(T-,v)dA 


Required  for  the  integration  is  the  recurrence  equation 


.Clx-.vl 


(rn  -  I  )a' 


(J  -  1 
(m  - 


y)  r  1 

n  .L  a”*-' 


The  recurrence  is  started  with  the  initial  integrals  in  ttie  equations 


^OO  r 

J  =  -  - 


-*(x-iyl  • 


•r  -  ly  J« 

-«tx-ty)  pi 


J'oc  j  p-atx-ty)  pi 

_  p-*(x^.yi  ^  -  df  =  -  Ei(-6(x  -  ly))  (74) 

a  *:  J-oo  < 

The  recurrence  is  cycled  in  descending  order  if  (5(x  --  ly)  satisfies  the  inequality 

!(5(x  -  ly);  g  17  (75) 

Otherwise  the  recurrence  is  cycled  in  ascending  order  with  an  initial  approximation 
for  m  =  64. 

For  the  series  expansions  the  arguments  have  Chebyshev  spacing  and  interpolations 
are  made  with  1 1  -point  central  Lagrange  interpolation  Preliminary  computations  have 
established  a  matrix  of  coefficients  such  that  the  coefficients  for  progressively  increasing 
powers  of  the  argument  are  obtained  with  the  product  of  an  array  of  values  of  the 
function  and  the  matrix  of  coefficients. 

If  X  satisfies  the  inequality 


7U  (/ 18\*  U  IB 

8i/  \V  3:  /  1/  X 


then  £  and  6  are  given  by  the  equations 


6 

£  =  — 

7 


s\  X  /  U  X 


and  the  last  interval  can  be  jettisoned  without  significant  error.  Otherwise  e  and  6 
are  given  by  the  equations 


in  which  case  the  Reynolds  number  is  less  than  36 

The  components  of  velocity  for  the  line  flux  are  computed  by  the  following  subroutine. 


*  ^»  */»^** 


t  K  »  *  *  r  i  >  •  ;  »' 


I 


FORTRAN  Fl  BROl'TINE  FOR  LINE  FLL'X  PROFILE 


The  free-stream  velocity  L’  is  given  in  argument  AO.  and  the  kinematic  viscosity  e  is 
given  in  argument  AN.  The  coordinates  x.  y  are  given  in  the  arguments  A\.  A>  Ascending 
and  descending  recurrence  relations  are  used  in  the  evaluation  of  Fourier  integrals 
for  a  line  flux  profile.  The  components  u.v  of  velocity  are  stored  in  functions  ,  Fv 

NONSLIP  BOI  NDARY 

Exploratory  computations  with  LMEXPr  have  indicated  the  nature  of  the  flow  from 
a  line  flux. 

If  X  and  y  are  decreased  to  zero,  the  contribution  to  integration  extends  to  large 
values  of  k'  where  A  approaches  the  approximation 


At  the  limit  of  small  x  and  y  the  velocity  is  given  by  the  equation 

V  =  -  i  ~ 


0  (79) 


(x-0,y-0)  (80) 


'•I'*’ 

HV 

is 

n 

S* 


I 


which  is  an  efflux  upstream  and  an  influx  downstream 

When  X  is  increased  to  infinity  downstream,  the  integration  with  respect  to  k  makes 
a  significant  contribution  to  the  integration  only  where  *■  -  0  and  A  approaches  the 
approximation 


Then  the  velocity  is  given  by  the  equation 


\  47^1^1 


f.  y  - 

r^x^r 


(k-O)  (81) 


(x  -  «=)  (82) 


This  is  a  Gaussian  distribution  of  velocity. 

In  a  continuous  distribution  of  line  fluxes  the  efflux  and  the  influx  at  the  leading 
edge  would  tend  to  cancel,  ft  is  assumed  that  for  integrated  distributions  the  Gaussian 
profile  is  adequate 

The  formula  for  velocity  does  not  meet  a  nonslip  boundary  condition  because  the 
velocity  at  y  =  0  varies  inversely  as  the  square  root  of  distance  downstream  However, 
a  continuous  distribution  of  line  fluxes  may  be  integrated  to  give  a  constant  velocity. 
As  a  consequence  of  the  equation 


r  .  r  -I  ’  «  1* 

I  ^  / - ^  \  x  -  n  ^ 

./ovavx-a  "o 


the  integrated  velocity  is  given  by  the  equation 


7T  Jo 


v(x  -  n ) 


(y  =  0)  (84) 


as  required  by  the  nonslip  boundary  condition  Tlius  the  strength  of  line  flux  per  unit 
displacement  is  given  downstream  by  the  expression 


(85) 


Ttie  Gaussian  distribution  is  equivalent  to  the  dispersion  of  a  line  pulse  of  vorticity 
which  has  diffused  outward  for  a  time  equal  to  the  expression 


T  -  a 
U 


(86) 


and  has  been  swept  downstream  with  the  speed  U . 

FOI  RIER  TRANSFORM 

Let  f(r)  be  a  periodic  function  cf  x  Then  it  can  be  approximated  by  the  Fourier 
series 


ioo V  a^cosnx  ^  ^  6„sin  nx 

n=l  n=I 

Inasmuch  as  cosines  are  even  and  sines  are  odd,  they  satisfy  the  equation 

■»  +  fr 


cos  kx  sin  mx  dx  =  0 


(87) 


(88) 


Application  of  the  Euler  theorem  to  the  following  integrals 


ij  dx  ±  I  J 


i{k  m)  _ 

*  1 

-n 

I 

3 

=  0  (fc  -  m)  (89) 


shows  that  the  trigonometric  functions  satisfy  the  orthogonality  relations 


/" 


cos  kx  cos  mx  dx 


sin  kx  sin  mx  dx  =  0 


Application  of  the  addition  theorem  leads  to  the  equations 


r 


dx  =  2tt 


J+»r  n*n 

cos^kx  dx  -  I  sin^mx  dx  =  n 

-  ir  t  *  -  w 

The  mean  square  error  for  the  Fourier  series  is  given  by  the  equation 

2  p  -t  IT  <M  OO 

=  —  i|<io  S  ci„cosnx-r  ^  6„sin  nx  - /(x)!^  dx 


(90) 

(91) 

(92) 

(93) 


Differentiation  with  respect  to  the  coefficients  shows  that  for  least  squares  error  the 
coefficients  are  given  by  the  equations 


/(x)  cos  nx  dx 


/(x)  sin  nx  dx 


(94) 


Because  of  the  symmetry  of  the  trigonometric  functions,  the  trigonometric  series  is 
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piven  by  the  equation 

n  =  +  oc 

f(r) 


«  n  = 

^  -L-  V" 

277  ^ 

n  =  -  »  ^  ' 


n  :  *  oc  /*■*  n 


f(s)  cos  ti{t  -  s)  ds  ~  —  '' 

2?:  ~ 


/(^■)  < 


in(x  -  s) 


ds 


n  -  -  OP  -n 


The  substitutions 


:t7x 


:7TS 


(95) 


(96) 


expand  the  range  of  approximation  from  27t  to  L.  The  functions  f{x)  and  f(s)  are 
replaced  by  the  functions  F(x)  and  /'(•s)  as  expressed  b>’  the  equation 


F(x) 


2  n=  +  oc. 


The  substitution 


^frn ,  V 

F(s)e'  L 


inn 


ds 


K  = 


t97) 


(98) 


replaces  summation  with  respect  to  n  by  integration  with  respect  to  v  in  the  limit  as 
L  -  «  The  approximation  of  summation  by  integration  requires  the  functions  to  have 
properties  of  integrability.  The  Fourier  series  is  approximated  by  the  Fourier  transform 


/’(x)e'‘"dx 


F{x)  =  j  A(k)  die 

«/  -eo 


(99) 


(100) 


In  a  multidimensional  space  there  is  a  Fourier  transform  for  each  coordinate.  Each 
transform  in  a  series  of  transforms  can  be  applied  to  the  amplitude  of  the  previous 
transform. 

The  two-dimensional  transform  is  given  by  the  equations 


^{q,  /S)  =  J*  J* F(x.  y)  fly 


F(x.  y)  =  J  J dex 

Let  yfx.y.  z)  be  defined  by  the  equation 


dp 


if(x.  y.  2)  =  >t{a.  P)  e'V“ 

Then  is  a  solution  of  Laplace's  equation 

dx^  dy^  dz^ 

wherever  2  ==  0  The  derivative  of  ^  with  respect  to  2  is  given  by  the  equation 
difi 

dz  ^ 


(101) 

(102) 

(103) 

(104) 

(105) 


where  the  sign  is  -  for  z  0  and  the  sign  is  *  for  2  '  0  The  difference  in  the  derivative 
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on  opposite  sides  of  the  plane  r  ^  0  is  -iTru  where  a  is  the  source  density  on  the  plane 
The  amplitude 


2  V  o'  -  A{cx.  fi) 


(106) 


is  the  amplitude  for  Tttct  For  a  unit  source  at  the  origin  ttie  amplitude  is  given  by 
the  equation 


A(n.p)  = 


1 


iir  \  a  -  p- 


(107) 


Two  solutions  of  Laplace's  equation  are  identical  to  within  an  additive  constant  if  they 
have  the  same  normal  derivative  on  a  boundary  Thus  the  inverse  of  distance  is  given 
by  the  equation 


1 


1 


1 


1 


\  X 


f  C//? 

\  O'  - 


(108) 


Application  of  the  Fourier  transform  to  the  real  exponential  factor  in  ttie  integrand 
leads  to  the  equation 


1 


1 

O-rr^ 


1 


^ \{ax  *0V  > 


(in  dp  dy 


(109) 


The  evaluation  of  this  Fourier  integral  can  be  completed  after  a  transformation  of 
coordinates 

Let  K.d.<p  be  polar  coordinates  with  polar  axis  in  the  direction  of  r  Then  a.p.y 
are  given  by  the  equations 


a  =  K  sin  6  cos  0  )?  =  x  0  sin  0 

da  dp  dy  -  ic*  sin  6  dk  dO  d0 

and  integration  with  respect  to  6.  0  leads  to  the  equation 


y  =  K cos  9  (110) 

(111) 


,  tJtrcostf  *  2 


Sin  0  dK  do  d(f> 


2  f  sin  KT 

-  I  -  dK 

T(  A  T 


(112) 


This  equation  confirms  the  existence  of  a  Fourier  transform  for  the  potential  of  a 
pole 

Let  r.  9.  0  be  polar  coordinates  with  polar  axis  in  the  direction  of  ic  Then  x.  y.  z  are 
given  by  the  equations 


X  =  r  sin  9  cos  0  y  =  r  sin  9  sm  0 

dx  dy  dz  -  sin  9  dr  dO  d(t> 

Another  solution  of  I.aplace's  equation  is  given  by  the  equation 


z  -  r cos  9 


d  /  I 


(') 


cos  9 


Its  Fourier  amplitude  is  given  by  the  equation 


XKrconB  «.2 


Sin  0  dr  dO  d(f> 


(113) 

(114) 

(115) 

(116) 


Integration  with  respect  to  9,  0  is  completed  with  an  integration  by  parts  to  give  the 
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t'quation 


_2_  r( 

Jo  i 


cos((cr)  - 


sin(A:r)  ^  dr 
icr  )  r 


(117) 


Integration  with  respect  to  r  is  completed  with  an  integration  by  parts  to  give  the 
equation 


if  sin(icr)  i 

2n^K  [  KT  Jq  2rT^ 


(118) 


This  equation  confirms  the  existence  of  a  Fourier  transform  for  the  potential  gradient 
of  a  pole. 

The  gradient  of  inverse  distance  is  given  by  the  equation 


,  1  \  ri  -  i/j  '  'k  if  r  r  ai  -  iSj  -*  vk  ,  .  , 

-  = - - 3  -  - r,  ^ ^  rf/S  dy 

^  2TT'  j  j  j  a--  7^ 


(119) 


The  velocity  of  a  point  vortex  of  unit  strength  is  given  by  the  equation 

v=-lk.v(-)=  j - - 3  =  -—1  f  f  (120) 

4n  \rl  4n  1^2  ^  J  J  J  a®  -  /S"  - 

This  equation  confirms  the  existence  of  a  Fourier  transform  for  the  velocity  of  a  point 
vortex 

The  circulation  around  a  circle  at  coordinates  r,  9  is  given  by  the  equation 


vdr  =  - 


,  sin®0 


(121) 


A  change  from  spherical  polar  coordinates  to  cylindrical  polar  coordinates  and 
integration  with  respect  to  2  with  the  aid  of  the  equation 


r“  dz  ^ 

J-®  \  *  z^  2 


(122) 


confirms  that  the  circulation  is  unity  around  a  line  vortex  of  unit  strength  and  of 
infinite  length 

INCOMPRFSSIBLE  FLITD 

Let  r  be  a  position  vector  in  physical  space  and  let  k  be  a  position  vector  in  wave 
number  space  The  position  vectors  are  defined  by  the  equations 


r  -  xi  yj  *  ck 


K  -  <»i  -  /?j  yk 


(12rt) 


where  i.j  k  are  orthoconal  unit  vectors 

Let  the  velocity  vector  v  at  the  position  vector  r  be  expressed  by  the  equation 


J  A(K)f’” 


where  k  is  a  vector  in  wave  number  space.  dK  is  a  volume  clement  in  wave  number 
space,  and  .\(k)  is  the  amplitude  of  (he  Fourier  component  with  wave  number  k  The 


♦4  '  *. 


*  .  I  t 


amplitude  is  given  by  the  equation 


A(k) 


1 

8n^ 


v(r)  c  "  dr. 


(125) 


where  ;dr  is  a  volume  element  in  physical  space  The  gradient  of  velocity  is  given  by 
the  equation 


Vv  =  i  J  kA(k)  c"  ■■  'ofK' 

the  divergence  of  velocity  is  given  by  the  equation 


V-v  =  i 


and  the  curl  of  velocity  is  given  by  the  equation 


JkA 

the 

i  j  k  A(k)  c  '*  '■ 


Vxv 


A(k)  c’*  ■■  dK 


The  Laplacian  of  velocity  is  given  by  the  equation 


VVv 


A{k)  c’*  ''  dK] 


The  amplitude  k  A(k)  of  the  divergence  V-v  is  given  directly  by  the  equation 

1 


k  A(k)  =  ^  V  vc-"^  dr\ 


(126) 

(127) 

(128) 

(129) 

(130) 


If  the  divergence  is  zero  for  every  r  in  physical  space,  the  amplitude  is  zero  for  every  k 
in  wave  number  space  The  amplitude  satisfies  the  orthogonality  equation 

kA(k)  =  0  (131) 


If  a  nonzero  divergence  did  happen  to  occur  in  an  incompressible  fluid,  a  pressure 
pulse  would  be  created,  and  the  nonzero  divergence  would  be  dispersed  by  a  potential 
flow 

Terms  which  are  quadratic  in  velocity  are  expressed  by  double  integrals  Let  k,  and  Kg 
be  variables  of  integration  in  the  double  integration  Incompressibility  is  expressed 
by  the  equations 

k,  A(k,)-0  (k,  -  K2)  A(k,  -  Kj)  =  0  K2  A(Ka)  =  0  (132) 


When  terms  with  the  same  wave  number  are  collected  in  the  integration  the  amplitudes 
are  collected  in  a  convolution  The  product  of  velocity  and  its  gradient  is  given  by  the 
equation 


v  Vv  ^  '  J'  J'  A(k2)k,  A(k,  -  Kg)  c"'  dK,i  c/KjI 
The  square  of  the  Gradient  is  given  by  the  ecjuation 

(Vv  Vv)  -  '  J  I  AIk^)  K,K,  A(k,  Kg)  dK,  dKgl 
This  scalar  is  ttie  divergence  of  the  vector 


(133) 


(134) 


I 


A(K2)- 


K,K,K 
Ki  K, 


'  A(k, 


Kg)  ^  dKillf/Kji 


(135) 
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The  gradient  of  pressure  is  given  by  the  equation 


it  yj  y« 

AvKa)-— ^ ^ ^ -A(k.  -  K,)  t**' idK, !  dK,|  (136) 

J  *^'iKi 

Substitution  in  the  equation  of  motion  leads  to  the  integro- differential  equation 

--  A(k,  )  -  -  I  fi-Ki A(k,  )  -  I  I  A(k,)  (  I - )k,  AfK,  -  Ka)  ic/Kal  -  —  K.  K, A(k,  )  (137) 

dt  J  \  K^  K,  /  p 

which  expresses  the  evolution  of  the  spectral  representation  of  the  velocity  in  an 
unbounded  fluid 

The  rate  of  change  of  amplitude  is  expressed  as  the  sum  of  three  terms  The  first 
term  replaces  jr  m  the  Fourier  transform  with  x-  I  t  The  second  term  expresses  the 
influenc-e  of  amplitude  at  other  points  in  wave  number  space.  The  third  term  gives 
the  rate  of  viscous  dissipation. 

The  rate  of  change  of  AtK,)  is  influenced  by  the  presence  of  A(>:a)  in  the  vicinity 
of  K,  In  the  integrand  tlie  postinultiplication  of  A(Ka)  by  a  tensor  eliminates  any 
component  of  AIk,)  in  the  direction  of  k,  Other  terms  outside  the  integral  are  the 

products  of  scalars  and  A(k,)  Thus  the  integro- differential  equation  preserves  the 

orthogonality  of  k,  and  A(>c,) 

There  is  no  contribution  to  the  integration  where  Kj  is  orthogonal  to  k,  and  AIkj) 
IS  parallel  to  k,  There  is  no  contribution  to  the  integration  where  Kj  is  collinear 
with  K,  and  A(k,  Kj)  is  orthogonal  to  k,.  If  Kj  is  on  the  perpendicular  bisector  of  k, 
then  interchange  of  Kg  and  k,  -  k,  leaves  amplitudes  the  same  but  reverses  the  sign. 
There  is  no  contribution  to  the  integration  by  integration  along  the  perpendicular 
bisector  of  k,  There  is  a  maximum  contribution  to  the  integration  when  Kj  is  on  the 
perpendicular  to  k,  through  the  lip  of  k, .  The  influence  of  AlKg)  is  a  pattern  which 
is  crossed  by  nodal  lines 


DlFFl'SION 

During  unbounded  evolution  the  velocity  deviates  gradually  from  the  boundary 
conditions  at  the  surface  of  a  plate  Velocity  is  injected  gradually  into  the  stream  to 
maintain  the  boundary  conditions  It  is  only  at  the  surface  of  the  plate  that  velocity 
IS  injec  ted  Fverywhere  else  the  evolution  of  velocity  is  free.  Let  the  plane  with  plate 
be  divided  into  an  equally  spaced  grid  Over  each  grid  point  there  is  a  sine  quotient 
function  Fach  sine  quotient  function  is  unity  at  its  own  grid  point  and  is  zero  on 
every  end  line  which  does  not  pass  through  the  grid  point  An  analytic  function  can 
be  expressed  as  a  series  in  sine  quotient  functions  The  coefficients  of  the  terms  in 
the  sc-nes  are  just  tfie  values  of  the  function  at  the  grid  points  The  terms  for  defect 
m  velocity  are  finite  only  for  grid  points  within  the  area  of  the  plate  The  sine  quotient 
function  at  any  grid  point  can  be  expressed  by  a  Fourier  transform  The  amplitude 
of  ttie  Fourier  transform  is  a  rectangle  of  constant  density  m  wave  number  space. 

l.cl  X,  y.  z  be  Cartesian  coordinates  relative  to  a  grid  jioint  with  x  positive  downstream, 
witti  y  positive  across  ttie  stream,  and  with  c  positive  perpendicular  to  the  stream 
Let  Q(a.fi)  be  the  Fourier  arrifilitiide  for  the  sine  quotient  function  for  the  grid  point 
and  hi  A  be  ttie  defect  of  velocity  at  ttie  grid  pcjint  The  velocity  is  injected  into  the 
laminar  sublayer  wtiere  Itie  aliscnce  of  convection  reduces  the  equation  of  motion  to 


The  deflect  iii  boundary  condition  is  concentrated  initially  against  the  plate  surface 
and  then  diffuses  outward  The  velocity  for  a  defect  A  is  given  by  the  equation 

V  -  A  "  c  “^du  (139) 

\  4l  t 

in  which  the  velocity  is  proportional  to  the  complementary  error  function.  That  this 
IS  a  solution  of  the  diffusion  equation  can  be  verified  by  direct  substitution  if  c  ^  0. 

The  error  function  is  s\nimetric  with  respect  to  the  plane  z~  0.  The  f'ourier 
amplitude  of  the  error  function  is  the  sum  of  an  integration  in  the  range  -<»  <  c  <  0 
and  an  inlearalion  in  ttie  range  0  v  c  The  integrals  in  ttie  two  ranges  are  complex 

coiijUi.'al  es  Inti'afat  ion  by  parts  m  cacti  range  and  cancellation  of  complex  conjugates 
h‘ad  to  the  equation 


du  dr 


*  y  \  -  7T  - 


+  V  ft 


du 
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The  ftmphtude  of  ttie  error  function  is  an  even  real  function  of  y  The  velocity  is  given 
by  ttie  equation 

•t>  Vf( 
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\  rr 


^  Q(o./S) 
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*■’  ■  ut 


ly 


da  d/S  dy 
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To  verify  this  equation  the  error  function  can  be  replaced  by  its  absolutely  convergent 
power  series,  then  term-by-term  integration  leads  to  the  series  for  the  arcsine  of 
unity  The  equation  gives  the  initial  velocity  for  the  defect  A  The  velocity  continues 
to  evolve  thereafter  in  accordance  with  the  full  equation  of  motion  for  unbounded  flow. 

DISCUSSION 

.•\n  analysis  of  the  complete  Navier- Stokes  equation  for  Reynolds  numbers  less  than 
1218  has  been  given  by  Scliwiderski  and  bugt’  They  found  overshoot  in  the  velocity 
profile  and  the  overshoot  increased  with  Reynolds  number  There  does  not  seem  to 
be  any  confirmation  for  the  overshoot.  There  is  no  overshoot  in  the  Blasius  profile. 
Validation  of  the  Blasius  theory  for  Reynolds  numbers  more  than  100000  is  given  by 
the  excellent  aereernent  between  measurement  and  theory  in  Figure  7.9  in  the  text 
by  Schlichtmg’  It  would  lie  nice  to  have  more  experimental  data  for  a  smaller  Reynolds 
number 

The  flux  of  momentum  across  any  cross  section  is  equal  to  the  drag  upstream  of 
the  cross  section  The  flux  of  mornentum  is  given  by  the  integral 


u  ( f ' -  u)  dy 


142) 


where  u  is  the-  instantaneous  velocity  in  the  fluid  The  instantaneous  velocity  can  be 
replaced  by  mean  vcdoc  ity  only  m  the  case  of  laminar  flow  Otherwise  the  mean  square 
of  the  instantaneous  velocity  is  greater  than  the  square  of  the  mean  velocity  by  the 
mean  scpiare  of  the  fluctuation  m  velocity 


Heretofore  it  has  been  assumed  that  the  blendine  function  for  a  turbulent  boundary 
layer  is  symmetric  witti  respect  to  the  middle  of  the  boundary  layer,  and  is  zero 
outside  the  thickness  of  the  boundary  layer  It  is  more  likely  that  the  blending  function 
diminishes  gradually  with  distance.  An  empirical  formulation  is  compared  with 
experimental  data  in  the  following  figure 


<5 


where  the  curve  represents  Equation  (35)  and  the  circles  are  from  a  tabulation  in  the 
report  by  Thompson®. 

The  assumption  that  the  horizontal  component  of  velocity  is  constant  everywhere 
in  the  free  stream  leads  to  a  vertical  component  which  persists  to  infinite  distance 
from  the  boundary.  However,  the  distribution  of  vertical  velocity  initiates  a  jet  and 
diffusion  of  vorticity  disperses  the  jet  in  the  free  stream. 

At  the  leading  edge  of  a  boundary  layer  the  flow  is  laminar,  but  the  flow  becomes 
unstable  w'here  the  Reynolds  number  is  520,  and  the  flow  becomes  fully  turbulent 
where  the  Reynolds  number  is  500000.  For  a  free-stream  flow  at  20  knots  the  laminar 
flow  IS  unstable  at  only  one  twentieth  of  a  millimeter  downstream,  and  the  flow  is 
fully  turbulent  at  five  centimeters  downstream.  Any  laminar  flow  in  the  boundary  layer 
of  a  full-scale  ship  is  insignificant. 

CONCLUSION 

It  is  concluded  that  the  most  practical  representation  of  the  mean  velocity  is  a 
generalization  of  the  empirical  formulation  by  Schlichting. 
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